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Abstract 

The sensitivity of the evolution of quantum uncertainties to the choice of the initial 
conditions is shown via a complex nonlinear Riccati equation leading to a reformulation of 
quantum dynamics. This sensitivity is demonstrated for systems with exact analytic solu¬ 
tions with the form of Gaussian wave packets. In particular, one-dimensional conservative 
systems with at most quadratic Hamiltonians are studied. 


1 Introduction 

The study of time-dependent (TD) phenomena is important in all areas of physics. In the 
framework of the time-dependent Schrodinger equation (TDSE) there are integrable model 
Hamiltonians which allow for analytical solutions. Examples related to the harmonic 
oscillator (HO) system are reviewed in [1]. In classical mechanics integrability, for a time 
independent Hamiltonian with n degrees of freedom, means that in phase space there are 
n independent constants of motion (see, e.g., [2,3]). 

Of course, for the conservative one-dimensional case the Hamiltonian system is always 
integrable. For TD Hamiltonians the situation is completely different because in this 
case the energy is not a conserved quantity and so even one-dimensional problems can be 
non-integrable [4], 

Wave packets (WPs) that propagate without dispersion were initially studied in quan¬ 
tum mechanics by Schrodinger [5]; these are packets of constant width that are usually 
discussed in contemporary textbooks on the matter and represent an early example of the 
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coherent states introduced by Glauber [6], Klauder [7], and Sudarshan [8]. That differ¬ 
ent initial conditions can lead to WPs with totally different properties of the width will 
naturally emerge from our investigation. This last point is significant as the HO is still 
essential in the high-precision measurements of systems with weak dispersion [9]. 

A different form of studying TD problems in quantum mechanics is via their propaga¬ 
tors, transforming an (initial) state at time t 0 into a (final) state at a later time t. Early 
on in quantum mechanics, Kennard [10] derived the propagators for the free particle, 
the motion of a charged particle in a uniform electric field, and the HO. This subject 
gained renewed interest when Feynman derived the propagators via his path integral for¬ 
mulation of quantum mechanics which has its origins in the classical principle of least 
action [11]. Since then an intensive search for propagators of Hamiltonian systems began 
(see also [12]). For (also multidimensional) Hamiltonians that are quadratic in posi¬ 
tion and momentum, together with those Hamiltonians associated with group-theoretical 
structures, there are established procedures to derive the propagators. As an example, 
the method of linear TD invariants introduced in [13,14] shall be mentioned because it is 
related to our approach for describing the time-dependence of quantum uncertainties, so 
that solving the problem is traced back to solutions of the classical equations of motion. In 
particular, the study of the parametric oscillator with TD frequency is very useful because 
of its applicability to quantum optics, plasma physics and gravitational waves [15,16]. 

In this work, we consider TDSEs with exact analytic Gaussian WP solutions. This 
is the first of a series of three papers that integrate our research. It is known that 
Gaussian WPs exist for any Hamiltonian that is at most bilinear in position and mo¬ 
mentum operators. In this case also exact analytic expressions for the time-evolution of 
quantum uncertainties exist; in particular for the HO, an oscillator with TD frequency 
uj(t) and the free motion. In the case of the parametric oscillator with u = oo(t), the 
Hamiltonian is no longer a constant of motion but the so-called Ermakov invariant still 
exists [17] 1 . Moreover, it is closely related to the Wigner function of the system [19], 
and also has been thoroughly used to study the parametric oscillator in classical and 
quantum physics [20-23]. In addition, the Ermakov invariant is relevant in studying the 
time-evolution of WPs associated to an oscillator with TD frequency. For example, this 
has been used in the determination of the oscillator’s squeezing properties when the fre¬ 
quency is a step function of time [24] . More recently, the Ermakov invariant has been used 
in a class of waveguide arrays where the refractive indices and second-neighbor couplings 
define the mass and the frequency of an analog parametric oscillator [25]. 

In section 2, the TDSEs for conservative systems with at most quadratic potential. 
Particularly the HO with time-independent frequency cuo are considered where the free 
motion can be obtained in the limit coq —> 0. The Gaussian WPs associated to these 
TDSEs are mainly defined by the coefficient of the quadratic term in the exponential, 
this coefficient is solution of a complex Riccati equation and includes information on 
the quantum uncertainties of position, momentum, and position-momentum correlation. 
Due to the quadratic term in the Riccati equation, one expects at least two different 

Where is an English translation of the work of Ermakov in [18]. 
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solutions (+ and — signs of the the square root). In addition, as the Riccati equation is a 
nonlinear (NL) differential equation it is sensitive to the initial conditions and can provide 
qualitatively quite different behaviour of its solutions depending on the choice of the initial 
conditions. However, these are related to different physical situations of the initial WP 
(e.g., minimum uncertainty WP or others). Because of these uncertainties also determine 
physical properties of the system, like tunnelling currents and ground states energies, 
these (measurable) properties are also affected by the choice of the initial conditions. 

In section 3 we show the form in which the initial uncertainties enter the solution of 
the complex Riccati equation. The relations between the Riccati equation, the Ermakov 
equation, the Ermakov invariant, and the aforementioned linear dynamical invariants are 
also established. There it is also shown how the potential (via the Hamiltonian), and thus 
the TD frequency of the parametric oscillator, can enter the treatment and consequently 
how our formalism can be extended easily to cases where no analytic solutions may exist 
for the classical equations of motion. A more detailed discussion will be provided in part 
III of our study concerning systems with explicitly TD Hamiltonians. 

The linearization of the Riccati equation leads to a complex Newtonian equation whose 
real and imaginary parts are not independent of each other and provide the TD parameters 
that essentially determine the propagator of the system. They also allow the rewriting 
of the Ermakov invariant in a form that the comparison with the Wigner function of the 
system is feasible. This is shown in sections 4, 5 and 6. 

In the final section, possible generalizations of our method are mentioned, especially 
the inclusion of dissipative effects which will be discussed in more detail in the part II of 
our study. There, we shall consider dissipative systems with friction that depends linearly 
on velocity. 


2 Riccati equations for conservative systems 


In the following we consider a one-dimensional oscillator with possibly TD frequency co{t). 
The corresponding TDSE is given by 


<9 T . \ 




possessing Gaussian WP solutions of the type 


^{x,t) = N(t ) exp < i 


y(t)x 2 + ~{p)x + K(t) 


( 1 ) 

( 2 ) 


where y(t ) = yR{t)+iyj(t ) is a TD complex function, x = x—(x)(t) = x—r]{t), rj(t ) denotes 
the position of the WP maximum following the classical trajectory and (p)(t) = mf](t) is 
the classical momentum. The purely TD normalizing factor N(t) and phase factor K(t ) 
are not relevant to the following discussion. The imaginary part of y{t ) is connected to 
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the WP width, or position uncertainty, via <r 2 (t) = ( x 2 ){t ) = yyyyy- Inserting WP (2) 
into the TDSE (1) one obtains, as equation for the WP maximum 2 , 


ij + co 2 (t)r] = 0 (3) 

and the dynamics of the WP width is defined by the complex Riccati equation 

C + C 2 + u 2 (t) = 0 (4) 

with C(t ) = — y(t ) = Cn(t) + iCj{t ) where an overdot denotes derivative with respect to 
time. 

The dynamics of the WP is completely determined by the solutions of the linear 
Newtonian equation (3) and the Riccati equation (4). 

To solve Eq. (4) we use the ansatz C{t) = C + V(t) where C denotes a particular 
solution of this equation that can be constant or TD. Thus, the inhomogeneous Riccati 
equation (4) turns into the homogeneous Bernoulli differential equation 

V + 2CV + V 2 = 0 . (5) 

It is well-known that this differential equation can be linearized via the transformation 
V(t) = k( t ) _1 to yield 

k - 2 Ck = 1. (6) 

The solution of this equation can be given in closed form and depends essentially on the 
particular solution C and the initial value K(to) = /to- For given C and initial condition 
k o, the solution of the Riccati equation (4) is fixed. 

To show the influence of the initial condition /to we consider, without losing generality, 
a HO with constant frequency uq. In this case, one immediately identifies two particular 
solutions to the Riccati equation (4), 

C± = ±iu 0 . 


The solution C + = iuio leads to the Gaussian WP solution with constant width, usually 
discussed in connection with the coherent states of the HO in quantum optics textbooks 3 . 
The second solution C_ = — iojq would replace the negative sign in the exponent of the 
Gaussian function by a positive one, leading to a function diverging for x —* ±oo, i.e., 
not to a normalizable WP. So, only one of the two possible mathematical solutions of the 
Riccati equation has physical relevance in this case. It should be noted that this situation 
can change when considering the complete solution C{t) = C_ +V(t). In this case, a 
positive imaginary part of V(t) could overcompensate G_ and thus allow for a second 

2 Depending on the particular time-dependence of u}(t), the solution of this equation might not even 
show an oscillatory or periodic behaviour. 

3 It was found already by Schrodinger [5] as a first example of what is now called a coherent state. The 
same WPs were retrieved by Glauber [6], Klauder [7], and Sudarshan [8] in their studies on the quantum 
coherence of light. 
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normalizable WP with TD width. Examples of this kind will be given in Part II of our 
research, on the topic of dissipative systems. 

For constant C±, the analytic solutions of the Bernoulli equation (5) take (with to = 0) 
the form 

1 p -‘2C±t 

V ± (t) — —— — -—----. (7) 

«±W K o + [1 - e~ 2C ^] 

As can be seen, the functions V±(t), and thus C(t), depend strongly on the initial condition 
k 0 . 



-4 -2 0 2 4 

C R 


Figure 1: Vectorial field of the Riccati equation for several initial conditions; the frequency is ujq = 1. 


The vectorial field of the Riccati equation for the HO with constant frequency is 
displayed in Fig. 1. It has two singular points at C± = ±fcao (the particular solutions) 
with periodic orbits around them. Due to the fact that the imaginary part of C(t) is 
inversely proportional to position uncertainty, all solutions in the lower half-plane are 
unphysical as they correspond to Cj(t) < 0. In addition to the above-mentioned choice 
of the particular solution, also the value of the initial condition k 0 can determine if a 
solution of the Riccati equation is physical or not. 


3 Riccati—Ermakov connections 

The connection of the initial condition k 0 to the initial properties of the WP (like its 
initial width) are demonstrated using a reformulation of the complex Riccati equation in 
terms of a real Ermakov equation. This can be achieved by introducing a new variable 
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a(t) via Ci(t) = ^k_. Inserting this into the complex Riccati equation (4) yields, from 
its imaginary part, the real part Cn(t) = ^j, leading to 


C\t) 


a(t) . 1 

a(t) a 2 {t) 


( 8 ) 


where a(t) is a real function of time and related to the WP width via a(t ) = 

From the real part of the Riccati equation one finally obtains 

a + u 2 (t)a = —. (9) 

cr 

Solving this Ermakov equation is equivalent to solving the Riccati equation. It has been 
shown by Ermakov [17,18] 1 that by eliminating uj (TD or time-independent) between 
equations (3) and (9) one obtains the dynamical invariant 4 5 * * * 


/ 




( 10 ) 


The quantum uncertainties of position, momentum and their correlation can be deter¬ 
mined directly by calculating the related mean values in terms of the Gaussian WP 
solution. They fulfil a closed system of coupled differential equations (see A), and can be 
expressed in terms of a{t) and a(t), or, in terms of the real and imaginary parts of C(t), 
respectively, as: 


<d(t) 


&xp (^ ) 


< i2 >« = = 2m ckr 


, o. . . mh 
(P 2 )(t) = — 


a 2 (t) + 


a 2 (t) 


mh C\{t) + Cj(t) 
"2 CAt) 


){[*.P]+)(<) = \{xp + pi)(i) = ^a(t)a(t) = 


( 11 ) 

( 12 ) 

(13) 


where Eq. (8) provides the relation between C(t) and a(t). Furthermore, it can be shown 
straightforwardly that also the Schrodinger-Robertson uncertainty relation 

4 The Ermakov equation and the corresponding invariant were reported by Ermakov in 1880 [17], over 
the years they have been reinvented several times by various authors [20,26,27], in a quantum mechanical 
context, particularly by [20]. However, the Ermakov system was actually discussed six years earlier in 
1874 by Steen [28]. But, since it was published in Danish in a journal not usually containing articles 
on mathematics, this work went unnoticed until the end of last century. An English translation of the 
original paper [29] and generalizations can be found in [30] 

5 Usually, this invariant is given without the factor Factor m came about in view of the fact that 

Ermakov considered the mathematical equation (3) whereas the corresponding Newtonian equation in a 

physical context must be multiplied by to. Since the resulting invariant has the dimension of an action, 

division by h provides a dimensionless constant. 
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is fulfilled. 


In order to obtain explicit expressions for the time-dependence of the uncertainties, 
the Ermakov equation (9) must essentially be solved for given initial conditions a (to) = «o 
and a(t 0 ) = do or the Riccati Eq. (4) for given kq. 

However, it is interesting to note that the solution of the Ermakov equation (9) can 
also be constructed knowing two linear independent solutions ?/i(t) and i^it) of the New¬ 
tonian equation (3). This can be achieved using the method of linear invariant operators, 
introduced by Manko et al. [13,16], and their relation with quadratic invariant operators, 
just as it is outlined in Appendix B. The solution of the Ermakov equation (9) can then 
be given in the form 

a(t) = J ^dg + + alnl{t) =F 2a 0 a 0 7] 1 (t)r]2(t ) 

I 2m 

= y ~Y + ^Mfa) T 2v xpo Vi(t)V2(t)] (14) 


where the two solutions of the Newton equation and their time-derivatives have the initial 
conditions 

Vi(to) = 0, r)i(t 0 ) = -, mito) = 1, mito) = 0- (15) 

m 

Given Eq. (14), the initial conditions of the Ermakov solutions can be expressed in terms 
of the position and momentum uncertainties together with their correlation function at 
initial time, i.e., 


ol o 



2 a. 


xpo 


hm a XQ 


This is consistent with Eqs. (11)-(13) for t = to: 


G^o) = <ri, a 2 Jt 0 ) = a. 


T 2 
PO’ 


a 


xp 


(to) = CT; 


XP 0 


There is also a relation between the initial condition Ko of the Riccati solution and the 
initial WP uncertainties. Using (11) and (13) one obtains 


C(t 0 ) = Co 


Q(n . 1 1 (J T7 ) n H 

i _—_+ i _ 

«0 «0 rn ° 2 xo 2ma l 


so that 

Vo = - = C 0 - C 0 . 

Note that all previous results are valid for TD as well as time-independent frequency to. 

As an example, the HO with constant frequency ujo is now considered. For the partic¬ 
ular solution C .|_ = iuio we have 


Uo = 


do 
Oi o 


+ i 



Kq = 





2 • 
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In terms of <7% 0 , cr^ and cr% po , the initial condition kq can be written as follows 



-4 ' -2 ' 0 ' 2 ' 4 


O' 

Figure 2: Vectorial field in phase space of the Ermakov equation for the parametric oscillator with 
constant frequency loq = 1. 


The vectorial field in phase space of the Ermakov equation for the HO with constant 
frequency u> 0 = 1 is displayed in Fig. 2. This has two singular points a = ±^=. The 
solutions describe orbits around these singular points; the ones on the right correspond 
to roots of a positive sign whereas the orbits on the left correspond to roots of a negative 
sign. 


Furthermore, as the integral curves in phase space are closed orbits, the functions 
a(t) and a(t) are periodic. Therefore, the uncertainties and the correlation function are 
also periodic. An exceptional case is obtained with the singular solutions for which one 
obtains 


o*{t) = 


h 


2 mcoi 


, ° 2 p (t) = 


Kmuj n 


&xp(t) 0 - 


These values correspond to the ground state of the HO and to the coherent state with 
constant width. Therefore, if the initial WP is not this coherent state, the uncertainties 
and the correlation function evolve periodically in time. 
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The analytic expression of the solution of the Ermakov equation is 


«(*) = \f 


? sin 2 u 0 t + al o cos 2 oj 0 t =F 


2<j x 


mojQ 


sin oj 0 t cos Loot 


x 0 


x 0 


= «o A / o + V 4 ) X 2 sill 2 u 0 t + cos 2 u 0 t =F ^ sin w 0 f cos w 0 f, 


"0 


cjq ao 


so that the uncertainties and the correlation functions are 
-2 

*■ I / I = - XIII (i lr\ I -t- 

#0 


°l(t) = 


So -2 J. , 2 2 J. , 2cr 3 -po 

Sill Wot + cS COS W 0 t + 


2,.,2 




o 


mujQ 


sin w 0 t cos w 0 t, 


) / \ ) ) ) ) ) ) 
cJt) = <j po cos uJot + m(jJ 0 a x sin — 2mcJo(T xpo sin uj 0 t cos(j 0 t, 


^xp(^) 




a. 


v 0 


mu>o<j. 


x 0 


sin 2 w 0 t + cr xpo cos 2 u 0 t. 


xpyj \2mtu 0 2 1 u ' xpo 

The corresponding results for the free motion are obtained easily using 

r sina; o t 

hm -= £, hm coscoqZ = 1 

COq 

and are presented in Appendix D. 


(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 


These general expressions lead to particular physical solutions once the initial values 
ct 2 q , <t 2 o and a xpo (or a 0 and d 0 ) are given. 


For example, consider a WP with parameter ao given in terms of constant frequency 
Wi = 0.5 as cio = This WP corresponds to the coherent state of a HO-Hamiltonian 
with frequency U\ (i.e., 1/a 2 7 ^ uo) ■ The initial uncertainties and correlation function 
are given, respectively, by 


h 


2 _ _ 

3:0 2 mwi ’ 


2 

So 


hmoui 


&xpo 0. 


The time-evolution of position and momentum uncertainties, with parameters m = 
h — 1, are displayed in Fig. 3a and Fig. 3b, respectively. Both uncertainties are periodic 
functions of time. In Fig. 3c, a 2 (f) and cr 2 (t) are plotted together with a straight line 
indicating the presence of squeezing of the WP (the uncertainty cr 2 (f) or cr 2 (f) takes a 
value below this line). In Fig. 3d the correlation coefficient, defined as 


Cor {t) 


| &xp (Z) I 
a x (t)a p (t) ’ 


( 22 ) 


is plotted. Note that Cor(t ) evolves periodically in time from zero to approximately 0.6. 


4 Riccati—Newton connections 


The Riccati equation (4) can be linearized replacing the complex variable C(t) by the 
logarithmic derivative of another complex variable, A (t) = A n(t) + according to 


C(t) 


a (ty 


(23) 
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crjd) 



cr?(t), t) 




Figure 3: Time-dependence of position (a) and momentum (b) uncertainties. In (c) both uncertainties 
are plotted together with a straight line in cyan colour determining if the WP is squeezed. In (d) the 
correlation coefficient Cor(t) is given as a function of time; in all cases we consider m = K = 1. 


Inserting this into the Riccati equation leads to the complex Newtonian equation 

A (t) +uj 2 (t)X(t) = 0. (24) 

Writing \(t) in polar coordinates as \(t) = the corresponding Riccati variable 

reads C(t ) = ^ This is in agreement with definition (8), providing 0(f) = 

holds, which looks similar to the conservation of angular momentum for the motion of a 
system with central forces. The validity of this conservation law can easily be shown by 
inserting the above form of C(t ) into the Riccati equation and looking at its imaginary 
part. It also follows from the Wronskian of the complex Newtonian equation, 

Xi(t)\ji{t) — = 1, (25) 

if real and imaginary parts of X(t) are written as A#(f) = a(t ) cos 0(f) and A j(t) = a(t ) sin 4>{t). 
Apparently, one of the components of X{t) is, up to a constant factor c, identical to the 
classical trajectory r]{t). For reasons that become clear later on, we choose A/(f) = cr]{t). 
Rewriting the Ermakov invariant (10) in terms of A i(t) and A i(t) instead of ?/(f) and 

r](t), it is obvious that the second term of the invariant reduces to = sin 2 0(f). 
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Therefore the first term of the Ermakov invariant must be cos 2 0(f), i.e., 

= (A/(t)a(t) - Aj(f)d(f)) = cos 2 0(f). 

This expression allows one to determine A^(t) in closed form up to a ± sign, 


A R (t) = ±ca 2 (f) ^(f) - ^y??(f)^ 


with c 



Using Eqs. (11) and (13), one can write in matrix form 

( Mt) \ _ ( Tcoi(t)a(t) ±^« 2 (f) A ( (x){t) A 

V Ar(f) ) V c 0 ){m) ) 

( T|^(f) ±|a 2 (f) A / <x)(f) A 

v c o j { < P )(t ); • 


( 26 ) 


(27) 


(28) 


Hence (x)(t) = r)(t) and (p)(f) = m?)(f) are still two linear-independent solutions of the 
linear second-order differential equation (3); but they are no longer orthogonal, whereas 
A n(t) and A/(f) are. In Eq. (14) two linear-independent solutions obeying Eq. (3) and 
the initial conditions (15) are needed to obtain the Ermakov variable a(t). One solution 
can always be the classical trajectory (x)(t) = rj(t) or A/(f) which differs from r/(f) only by 
a constant factor. As a second solution A n(t) can obviously be chosen. In cases where do 
or = Cr 0 vanish, also 77 (f) can be used because, according to (27), A#(t) and 77 (f) then 
fulfil (up to a constant factor) the same initial conditions. In the cases of the free motion 
(V(x,t) = 0) and the HO, with Gaussian WP solutions (also with oscillating width) one 
actually finds this situation, — = 0. 

J 7 QO 

Rewriting the Ermakov invariant (10) in terms of A#(f) and A/(f) as 


m 1 
2h~^ 



(29) 


and using (28), it can be expressed in terms of bilinear combinations of the mean value 
of position (x)(t) = r/(t) and momentum (p)(t) = mf}(t) multiplied by the corresponding 
conjugate uncertainties as 

^ [a 2 (f)(a;) 2 (t) - 2 a xp (t)(x)(t)(p)(t) + a 2 (t)(p) 2 (t)] . (30) 

Finally, writing a 2 (f) as a 2 = A 2 R (t) + A 2 (f) and using (28), one obtains an expression that 
is in agreement with Eq. (14). 
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5 Feynman kernel/propagator 


A different way of describing the time-evolution of a quantum system is to apply a prop¬ 
agator G(x,x']t,t') to an initial state T( x',t') that transforms this state into \l/(x,£) at 
position x and time t according to 


T(x, t) 


dx'G(x, x'] t, ^ , )\^ r (x , , t! = 0). 


(31) 


The integral kernel G(x, x';t, t') is also called Feynman kernel and was determined by 
Feynman using his path integral method [11]. 

In the cases discussed in this work, the propagator can be expressed easily in terms of 
real and imaginary parts of the complex variable A (t) that fulfils the complex Newtonian 
equation (24). Particularly for our case with Gaussian WP solutions of the TDSE, it can 
be shown that for an initial Gaussian state that can be written as 6 


T(x',F = 0) 



1/4 

exp 




, 0 P o 
+ 2 —x 
m 



(32) 


with po = (p) (to = 0), also the propagator has the form of a Gaussian function. The 
exponent of this function then obviously contains terms proportional to x 2 , x' 2 and xx' 
with appropriate TD coefficients. These coefficients can be easily expressed in terms of 
A n(t) and A i(t) and the propagator can be written in the form 


G(x, x'\ t, t' = 0) 


/ m A 

\2iriha 0 \i(t)) ^ 


irn 

2h 


Ar(t) 2 
A i(t) 




a R (t) 

A i(t) \ cxq / J 


(33) 


As the initial WP T(x / ,f / = 0) does not depend on x and t, the propagator G(x,x']t,t') 
itself must also fulfil the TDSE. Inserting G(x,x']t,t') into this equation leads to terms 
proportional to x 2 , x and independent of x. The coefficient of the latter ones provide the 
conservation law (25); the coefficient of the terms proportional to x cancel each other. 

From the ones proportional to x 2 , one obtains a Riccati equation for as variable 

(instead of C(t) = yyyy) that can be linearized to the Newtonian equation for A j(t) (which 
is proportional to rj(t) = (x)(t)). 

Inserting the propagator (33) into Eq. (31) and performing the integration now leads 

6 If the initial change of the WP width, and thus do, is different from zero, the term i (yy) in Eq. 

(32) has to be replaced by + i \ j x' 2 = C 0 x' 2 . However, in the examples discussed in this paper, 
this is not the case. 
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to a Gaussian WP in terms of x and A/(f), \n(t) instead of x and (x)(t), ( p){t ), i.e., 


^ wp(x,t ) = 


/ m \ V 4 

\ 7 rfi,/ V A^(f) + 


1/2 


exp 


ini 

2h 


^i(t) 

A/(f) 


a : 


a : 


*A,(f)) s 


A/(t)(Afi(f) + 


( 34 ) 


Comparing this with the WP in the form (2) shows that ^A/(f) = ?/(f) justifies our choice 
made above. With the help of conservation law (25), it is also easy to prove that the 
coefficient of the terms quadratic in x simply equals iy(t ) = (^?/(f)) = ffC'(i). 


6 Wigner function 

The Wigner function of the systems under consideration can be obtained from the WP 
solution (2) via the transformation 

W(x,p; t) = T. J°° dq ** ( x + i) >p ( x - |, i) e hr (35) 

and can be written as 

^ 1 [ 0 /"2/f 4” y“ii \ ~2 P 2 , 2?/R~~l , oa ^ 

H/n. P ;«) = ^exp|-2^——jx + (36) 

Using (8), (11), (12) and (13), this finally acquires the form 

W(x,p;t) = ^ ex P - 2a xp (t)x(t)p{t) + (J 2 x {t)p 2 {t)] \ . (37) 

Comparing this with the Ermakov invariant (30), one realizes that the Wigner function 
can be expressed as 

W(x,p-t) = -L e -^,p-,t) 
nh 

If we consider (x)(t) = 0 and (p)(t) = 0, this is reduced to 

W{x,p-,t) = (38) 

Tin 

Thus, when the expectation values of position and momentum are zero for any time, the 
Ermakov invariant I(x,p\ t ) explicitly defines the analytic form of the Wigner distribution 
W(x,p;t). Remarkably, in phase space the invariant I(x,p-,t ) describes ellipses that are 
rotating and thus their major and minor axes are changing in time. This behaviour is 
inherited by the maximum of the Wigner distribution which will describe also elliptical 
trajectories with the rotating axis directed along its maximum value. 
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Figure 4: The Wigner function for the parametric oscillator with parameters m = h= ojq = 1 where 
the initial state is a coherent state with parameter depending on the frequency ui\ = 0.5. This implies 
that cr^ 0 = 1 and a^ o = 2. We also use (x)(to = 0) = 1 and {p)(to = 0) = 2. The times used are 
t = {0, 7t/4, 7t/2, 37t/4, 57t/4, 37t/2, 77t/4, 27t}. 


In the most general situation, i.e. (x)(t) ^ 0 and (p)(t) ^ 0, the Wigner function 
displays the same behaviour; however, the maximum of the distribution will move in 
phase space according to the classical trajectory ((#)(£), ( p)(t )). 

As an example, we consider the Wigner function of an initial WP with its parameter 
depending on the frequency U\ and that is associated with a HO with constant frequency 
uj\. Thus, the initial uncertainties and correlation function are given by a 2 g = 2 mu \ , 
a po = WT 1 an< ^ ax Po = 0; respectively. The related Wigner function is plotted in Fig. 4 
for different times where we have used the parameters m = h = ouq = 1, or = 0.5 and 
the initial expectation values (x)(t 0 = 0) = 1 and (p)(to = 0) = 2. Notice that the 
maximum of the Wigner function describes an elliptic trajectory that is defined by the 
classical energy of the system 

= E d . (39) 

Furthermore, the probability distribution functions in position and momentum represen¬ 
tation can be obtained via 


T*(x, t)ty(x, t) 


|T(x,t )| 2 

\*M \ 2 


' —oo 
rco 


dp W(x,p;t ), 
dx W(q,p ; t). 


(40) 

(41) 


For the case under consideration, the probability distribution functions in position and 
momentum representations are plotted in Fig. 5. Note that the maximum of the distribu- 
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(a) 

Figure 5: (a) Probability distribution in position representation, (b) Probability distribution in mo¬ 
mentum representation. 


tion follows the classical trajectory and its width changes periodically in time according 
to the values of a x (t) and a p (t). 

Similar to the probability distributions t)T(x, t) = |^/(x,t)| 2 = p(x,t), the 

Wigner function also fulfils a continuity equation but now, as in the classical case, in 
phase space, i.e., 

^W(x, F ,t) + r -^W(x, P] t) - 9V ( X,t) %-W(x, F ,t) = 0 (42) 

at m ox ox op 

with p = p + (p)(t), x = x + (x)(t). 

Inserting (37) into this equation together with the TDHO potential one has at most 
bilinear or quadratic expressions in terms of x and p. The terms proportional to p 2 
yield Eq. (A.l) in Appendix A; those proportional to x 2 lead to Eq. (A.2) and those 
proportional to xp to Eq. (A.3). So, the set of three coupled first-order differential 
equations connecting the time-evolution of the uncertainties are regained in this way. 


7 Conclusions and perspectives 

It has been shown that the information about the dynamics of the quantum systems 
under consideration is not only obtainable from the TDSE, but also from a complex 
Riccati equation. Solution of this equation directly provides the information about the 
time-dependence of the quantum mechanical uncertainties of position and momentum 
and their correlation. As the time-evolution of these quantities essentially determines 
observable quantum effects like tunnelling currents (see Appendix C), the laws describing 
this evolution and the dependence on initial conditions are of major importance. The 
way the initial conditions enter these laws has been shown explicitly in Sections 2 and 
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3. Using the standard textbook example of the HO, it has been shown how a simple 
change of the initial conditions (WP with a width of the ground state Gaussian function 
vs. a width differing from this) can lead from a static situation (WP with constant width) 
to a dynamic one (WP with oscillating width). This is due to the fact that the static 
width refers only to a particular solution of the corresponding Riccati equation, whereas 
the oscillating width is related to the general solution. A possible second solution of 
the Riccati equation is discarded in the static case because it would lead to a diverging 
WP. In the case of the general solution, also this second particular solution can lead to 
physically-reasonable results (examples will be given in Part II which considers dissipative 
systems). 

Furthermore, the connections between the classical degrees of freedom and the corre¬ 
sponding quantum uncertainties have been shown in detail. Particularly from the linear 
invariants (see Appendix B), linearizing the Riccati equation (4) to the complex New¬ 
tonian equation (24) together with the conservation law (25) can show how quantum 
uncertainties can be expressed in terms of classical position and momentum and vice 
versa. 

The real and imaginary parts of the complex Newtonian equation also allow for the 
definition of the propagator or Feynman kernel of the system as shown in (33). In polar 
coordinates, the absolute value of A (t) is identical to the variable a{t) that transforms 
the complex Riccati equation into the real Ermakov equation and is directly proportional 
to the WP width. With the equation of motion for a(t) and the classical trajectory 
rj(t), it is possible to find a dynamical invariant, the Ermakov invariant (10), that is still 
an invariant in cases where the Hamiltonian of the system is not (e.g., for an oscillator 
with TD frequency oj(t) or for certain dissipative systems as will be shown in Part III). 
Moreover, this invariant allows one to immediately write down the corresponding Wigner 
function of the system. 

These results, obtained from standard quantum systems with exact analytical solu¬ 
tions, can however serve as the basis for generalizations. In Part II, it will be shown 
how the formalism can be extended to dissipative systems with linear velocity dependent 
friction forces. This will still be possible in the framework of analytical solutions. The 
next challenge will be systems (without and with dissipation) where no analytic solutions 
of the corresponding classical problem exist, i.e., Newtonian equations with potentials 
V(x, t) that cannot be solved in a closed form. In these cases, it is reasonably assumed 
that Newtons equation still applies. 

The same will be assumed for the Ermakov equation (corresponding to a complex 
Newtonian equation). As we know the connections between the classical variables and the 
corresponding quantum mechanical uncertainties (particularly the position uncertainty 
related to a(t)), from the problems with analytical solutions discussed in this paper, 
it should be sufficient to calculate numerically the solutions for rj(t) and fj(t) in these 
cases to derive the corresponding quantum uncertainties and properties like tunnelling 
currents and ground state energies. Besides, once r}(t), fj(t), a(t ) and a(t) are available, 
also A(£) = A n(t) + can be ascertained and, thus, the parameters defining the 
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propagator of the system and, via the Ermakov invariant (that might not be an invariant 
in these cases), the corresponding Wigner function determined. 

That the values of quantum mechanical properties like tunnelling currents are sensitive 
to the choice of initial conditions (and may change qualitatively according to this choice) 
can already be shown in analytic form considering dissipative systems with linear velocity 
dependent friction forces. This is dealt with in Part II. 


A Equations of motion of the quantum uncertainties 


For the parametric oscillator (for Aca(f) = 0 this also applies for the HO), the uncertainties 
of position and momentum and their correlation satisfy the following closed set of first 
order differential equations 


dt 

dt 

dcr xp (t) 

dt 


&xp(t), 


m 


= —2mui 2 (t)a xp (t), 




m 


muj 2 {t)a 2 x [t). 


(A.l) 

(A.2) 

(A.3) 


This system has an invariant given by Isr = cr 2 (t)a 2 (t) — cr 2 (t), corresponding to the 
Schrodinger-Robertson uncertainty product [32], 

To solve the system defined by (A.l), (A.2) and (A.3), it can be written as a third-order 
ordinary differential equation for cr 2 (t), 


d 3 v 2 x (t) , , 2 f4 .sdo%(t) , . / s du(t) 2 


+ 4 u 2 (t)- 


+ 4 u(t)——a x (t) = 0. 


(A.4) 


dt 3 ' y ' dt ' y ' dt 

This equation is known as the normal form of maximal symmetry because all linear or 
linearizable third-order ordinary differential equation can be transformed into (A.4), [33]. 

Inserting the expression 

h 

(A.5) 


^ 

into (A.4) shows that (A.5) is only a solution of (A.4) if a(t) is a solution of the Ermakov 
equation (9). From (A.5) and (A.l) the position-momentum uncertainty correlation a xp (t ) 
can be obtained immediately as 


o xp (t) = ^a(t)a(t) 


and from cr xp (t) and (A.3) it follows that 


<#) = '-f ( a 2 «) + 


(A.6) 


(A. 7) 
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B Ermakov solutions and linear invariant operators 


According to classical Hamiltonian mechanics, the integration of equations of motion 
becomes trivial if one can find a canonical transformation to variables that are constants. 
The same applies to the corresponding quantum mechanical operators. The linear TD 
operators X and P for the parametric oscillator are defined by the transformation 


( X \ _ ( 9i{t) hit) \ ( % \ 

[p J { 92{t) hit) ){p ) 


(B.l) 


which satishes the commutation relations of position and momentum operators, with the 
initial condition X(t 0 ) = x, Pft o) = p. The time-invariance implies that the functions 
gift) and fift) satisfy the coupled system of differential equations 

gift) = 

hit) = — -gift), 

m 

i — 1,2, which is equivalent to the second-order differential equation 

f,+u 2 (t)f, = 0. (B.4) 


(B.2) 

(B.3) 


With these invariant operators the most general quadratic invariant operator can be 
proposed as 


i 


AX 2 + BP 2 + C (xP + Px) 


(B.5) 


where A, B and C are constants still to be determined. As we already know a quadratic 
invariant for the system, the Ermakov invariant we assume that the invariant (B.5) is 
proportional to the operator corresponding to the Ermakov invariant, 7 , 


* 771 

” 2 h 


p 


aft )- aft)x + 


m 


x 


aft ) 


(B.6) 


Taking Eq. (B.5) to be equal to (B.6) by expressing X and P according to (B.l) in terms 
of x , p and comparing the coefficients of x, p and xp + px, one obtains for aft) 


aft) = ±Vmh \Aflft) + Bflft) T 2\C\f l (t)h(t )\ 1/2 , 


with 


\C\ = ±\lAB -- 2 . 


(B.7) 


(B.8) 


Using the initial conditions for f\ (t) and f 2 (t) and requiring that a 0 = a(t Q ) a 0 = d(t 0 ), 
one obtains 


m 


h 


l 


(y.r 


CXr 


A =xW + - 2 l b = -x, \c\ = 


mh 


ot o a o 
h 


(B.9) 


7 The usual definition of the invariant is without the factor . The conventional expression is obtained 
for m = h = 1. 
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The coefficients A, B and \C\ can be expressed in terms of the initial uncertainties and 
correlation function as 

A = fi2 a P0i B = Jp a x 0 i |C| = ff a xpo- (B.10) 

Obviously there is an ambiguity in sign in the expression for a{t). Considering the 
negative sign, the solution a{t) admits that a(t 0 ) and d(f 0 ) ma Y be different from zero. 
Considering the positive sign, a(t o) is still different from zero, but a(t 0 ) = 0. Therefore, 
in the more general first case, a(t) takes the form 

I 2m 

a(t) = \J— [crlJl{t) + a 2 x J 2 (t) T 2cr xopo f 1 (t)f 2 (t)]. (B.ll) 


C Quantum uncertainties and tunnelling currents 


The continuity equation describing the time-evolution of the probability distribution p = 
T*(a;, t) = |T(x, t)\ 2 in position space which corresponds to the TDSE (1) is given 

by 

d d 

dt P ^ Xl ^ + dx ^ Xl ^ = °’ (C- 1 ) 

where the velocity field v_(x,t) is defined as 


V-(x,t ) 


Jh_ 

2 im T(a;,f) 


JR’OM) 


h d T(x, f) 

-In- 

2 imdx T*(x, t) 


(C.2) 


and the corresponding current is p(x,t)v-(x,t)- 

The velocity V-(x, t) obviously only depends on the phase of the wave function T(x, t) 
and, therefore in our case of Gaussian WPs with complex coefficient C(t), also on Cn(t) = 
For these WPs the velocity field can be written explicitly as 

v- (x, t ) = pit) + ^jp-x (C.3) 

a(t) 

which shows that, apart from the classical velocity r)(t), also the tunnelling contribution 
depending on j s involved. Different initial conditions of the Riccati or Ermakov 
equation can lead to different time-evolution of a(t ) and thus to different tunnelling 
currents. 


D Uncertainties and correlations for the free motion 


In the case of the free motion (V = 0) the results corresponding to Eqs. (17) and (18)-(20) 
can be obtained easily using 


sin w 0 1 . 

Inn - = t, 

^ 0—>0 COq 


lim coswo t 

LJQ—tO 


1. 


(D.l) 
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The solution of the Ermakov equation can be written as 



+ < =F 2 a -^t 
m z m 


i±-d + 4 


«o / a 0 

and the uncertainties and correlation functions take the form 


(D.2) 


o*{t) 

m z m 

(D.3) 

°l(t) 

2 

(D.4) 


a 2 



= -5 °t + a xpo . 
m 

(D.5) 


For do = 0, or a xpo = 0, the well-known textbook results are regained. If these initial 

values do not vanish, two additional solutions for a 2 {t) are possible. However, one must 

carefully check if both mathematical solutions are also physically reasonable. For example, 

calculating a xp (t ) according to (A.l) by taking the derivative of (J 2 (t) would lead to 
2 

o xp {t) = —£ =)= a xpo . For t = 0, this would also allow for a xp (t 0 ) = —cr xpo which is 
obviously incorrect if a xpo 0. So, only the plus-sign provides a physically reasonable 
solution corresponding to a physical system with non-vanishing initial correlations (do 7 ^ 
or a xpo ± 0 ). 
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